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Abstract. We report on spectral and timing analysis of BeppoSAX data of the 13.6 s period transient X-ray 
[ pulsar EXO 053109-6609.2. Observations were carried out in March 1997 and October 1998, catching the source 

during a high and a low emission state, respectively. Correspondingly, the X-ray luminosity is found at a level 
Ijr^ • of 4.2 x 10 37 ergs _1 and 1.5 x 10 36 ergs _1 in the two states. In the high state the X-ray emission in the energy 

range 1-100 keV is well fitted by an absorbed power-law with photon index V ~ 1.7 plus a blackbody component 
with a characteristic temperature of ~ 3.5 keV. Moreover, we find an evidence of an iron emission at ~6.8keV, 
typical feature in this class of sources but never revealed before in the EXO 053109-6609.2 spectrum. In the low 
state an absorbed power-law with V ~ 0.4 is sufficient to fit the 1-10 keV data. During BeppoSAK observations 
EXO 053109-6609.2 display variations of the pulse profile with the X-ray flux: it showed single peaked and double 
peaked profiles in the low and high state, respectively. Based on these two observations we infer a spin-up period 
derivative of -(1.14±0.08) x 10" 10 ss' 1 . By comparing these with other period measurements reported in literature 
|' we find an alternating spin-up and spin-down behaviour that correlates well with the X-ray luminosity. 
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^ ' 1. Introduction and relative velocities along the neutron star orbit, which 

in most cases may results in regular X-ray outbursts near 



EXO 053109-6609.2 is a High Mass X-ray Binary (HMXB) 



the periastron (Type I outbursts). Alternatively, aperiodic 



' hosting a neutron star and a B-emission spectral type , , . . . , , , 

3 . . ° , TT . , , . _,. , „ outbursts occur which oltcn last longer than the neutron 

LU , (Be star Haberl, Dennerl, & Pietsch 1995; McGowan , .. , , , , , ... 

I ' , „«„\ ', . . , , star orbit. Inese are probably caused by matter ejection 

& Charles 2002). lire system is m the Large Magellanic ,„ . £ , . , , r , /rT , 

, , T ,,„ N '„, r -rms-,^* , i. , outflowing from the equatorial plane ol the Be star (lype 

Cloud LMC), ~ 17 away from LMC X-4, another hrgh n outbursts Stell white & Rosner 1986 Motch et aL 



mass X-ray binary pulsar (La Barbera et al. 2001). Be 



1991). 



stars are characterised by high rotational velocities (up nrmnn ^nnn n v ] • ,„ M , 

„ nM , „ . , , , ni . . . . , r EXO 053109-6609.2 was discovered m 1983 deep 

to 70% ol their break-up velocity), and by episodes ol ^v^oArn e n thi^i -v * ■ /n i n ±. 

, . , , n a, EXOSAT exposures of the LMC X-4 region (Pakull et 

equatorial mass loss which produce a temporary decre- n -,„„- , , „ , , „ -,^o^\ 1 , 

. ..... , , , , . , , . „ al. 1985: Pietsch, Dennerl. and Rosso 1989). lire lumi- 

tion disk around the star (Slettebak 1987, Okazaki & „ lri a fi /n , . Tr . .. 

. „„„ 1N m ,. . . ... nosity was ~ 6 x 10 JD ergs (0.15-4 keV: assuming a dis- 

Megueruela 2001). Ihese X-ray binaries m which a neu- 



tron star is orbiting in a relatively wide orbit with mod- 
erate eccentricity are the most numerous among HMXBs. 
In the LMC, more than half of the confirmed HMXBs are 
variable or transient sources consisting of a neutron star 
with a Be companion (Haberl, Dennerl, and Pietsch 2003; 
Sasaki, Haberl & Pietsch 2000; Haberl & Sasaki 2000). 



tance of 50 kpc). Subsequently, the coded mask X-ray tele- 
scope SL2 XRT flown on board of the shuttle Challenger 
between July and August 1985 detected a second out- 
burst (Hanson et al. 1989), with a source luminosity of 
- 1 x lO^ergs- 1 (2-10keV). 

The source was then monitored from June 1990 to 
July 1994 with the ROSAT PSPC (Haberl, Dennerl, & 



Variable X-ray emission is often observed which likely _. , . .„„ rN ... . . ... . , T .„q 7 , 

,, » ,. , . ,. . J .. „ . , . Pietsch 1995). Alternating high states (L x ~ 10 J 'ergs i ) 
results from the large variations m the Be wind density , . , , T ■•^■w _n . v n 
Z_ and low states [L x ~ 10 ergs ) m the X-ray flux were 
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from March to May 1993 with an average luminosity 
- 2.4 x lO^ergs" 1 (0.1-2.4 keV). 

Dcnnerl et al. (1996) reported on a ROSAT obser- 
vation that took place in October 1991 and led to the 
detection of a spin period of 13.67133(5) s with a period 
derivative of (1.5±0.1) x 10~ 8 ss _1 (calculated during their 
observation). Haberl et al. 1995 proposed an orbital pe- 
riod of about 600-700 days and an orbital eccentricity of 
e ~ 0.4 — 0.5. Under the assumption that period changes 
are caused by Doppler shifts, Dennerl et al. (1996) cor- 
rected the proposed orbital period, finding an orbital so- 
lution with P or t, = 25.4 days and an eccentricity of e ~ 0.1. 

Timing analysis of the BeppoSAX data obtained in 
March 1997 when the system was in high-state, revealed 
coherent pulsations at a period of 13.67590(8)s and two 
different period derivatives: a short-term period deriva- 
tive, calculated during the 2 days of BeppoSAX observa- 
tion, of Pi oc = (3.7 ± 0.5) x 10 _9 ss _1 , and a secular pe- 
riod derivative, calculated from a comparison with a pre- 
vious ROSAT measure of P sec = (3.67±0.05) x 10~ n ss' 1 
(Burderi et al. 1998). 

Last observation of EXO 053109-6609.2 was a deep 
XMM-Newton observation of LMC field of October 2000 
(Haberl, Dennerl and Pictsch 2003) in which EXO 053109- 
6609.2 was found to be the brightest source in the field 
(7 x 10~ 12 ergcm~ 2 s _1 in 0.2-10 keV range, corresponding 
to a luminosity of about 2.1 x 10 36 ergs _1 ). 

Across different observations this source varied its in- 
tensity up to a factor of about 10. Nevertheless, notwith- 
standing most of Be X-ray binary systems show a quies- 
cent state, this source has not yet been detected in quies- 
cence but only in a low-luminosity state with a minimum 
luminosity of ~ 10 35 - 10 36 ergs _1 . 

In this paper we report on spectral, spin period and 
pulse profile changes with the X-ray luminosity. We finally 
compare our results with previous findings. 
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Fig. 1. Epoch folding search of the data from the low state 
observation of 1998. 



2. Observations 

BeppoSAX observed EXO 053109-6609.2 with its 
narrow-field instruments (NFI; Boella et al. 1997a) from 
10:48:16 of 1997 March 13 to 08:17:50 of 1997 March 15 
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Time (TJD) 

Fig. 2. Evidence of a quadratic component in the phase- 
fitting of the high-state observation (linear component, 
corrisponding to the local P was removed in order to bet- 
ter show the quadratic residuals). 

and from 22:39:47 of 1998 October 20 to 08:05:31 of 1998 
October 22. 

Spectral coverage of BeppoSAX extends over a wide 
energy range using four different instruments. In our anal- 
ysis we use just two coaligned instruments: Medium- 
Energy Concentrator Spectrometer (MECS) consisting of 
three units operating in the 1-10 keV range (Boella et al. 
1997b) and the Phoswich Detector System (PDS) with 
four scintillation units operating in the range 10-200 keV 
(Frontera et a. 1997). The High Pressure Gas Scintillator 
Proportional Counter (HPGSPC, 7-40 keV; Manzo et al. 
1997) was not available; and for what concerns the Low 
-Energy Concentrator Spectrometer (LECS, 0.1-4 keV; 
Parmar et al. 1997), in the low state observation there 
are not enough photons to be used, while for the high 
state the off-axis response matrix is not accurate enough 
(actually LECS data were not fully consistent with the 
MECS data in the overlapping range ~l-5keV). Both in- 
struments used have a field of view less than one degree 
but only MECS has imaging capabilities. 

We extracted EXO 053109-6609.2 data from 4' radius 
circular regions in the MECS field of view, centred on 
the maximum of the point spread function of the source. 
The data were background subtracted by extracting in 
the same field of view background photons from a similar 
circular region of 4' radius, centred far enough from any 
sources present in the image and at the same off-axis po- 
sition of EXO 053109-6609.2; the arrival times were then 
corrected for the barycentre of the solar system. Owing 
to the problems descrived above, we used only MECS and 
PDS data for the first observation: we could use the PDS 
data because the contaminating source LMC X-4 was in 
a low state during that observation (Burderi et al. 1998). 
On the other hand, during the second observation LMC 
X-4 was in a high state (Naik & Paul 2003). In the PDS 
range, the relative fluxes between EXO 053109-6609.2 and 
LMC X-4 are: Fexo/Flmc ~ 9 for the first observa- 
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Fig. 3. First and second panels: pulse profile of March 
1997 BeppoSAX observation of EXO 053109-6609.2 in 
PDS and MECS bands aligned in phase, respectively. 
Third panel: pulse profile of October 1998 BeppoSAX 
observation of EXO 053109-6609.2 in MECS band, while 
the source was in low-state (not phase aligned with pre- 
vious observation) . 



tion, and Fexo/Flmc ~ 0-2 for the second (latter value 
was inferred extrapolating EXO 053109-6609.2 low state 



spectrum in the PDS range). Therefore, because during 
the second observation LMC X-4 was greatly more lumi- 
nous than EXO 053109-6609.2, and given that PDS does 
not have imaging capabilities, we only used MECS data 
in this case. 

MECS spectra, accumulated from the same circular 
regions used for the event files, were re-binned with the 
BeppoSAX grouping file: mecs_5. grouping; PDS spectra 
were re-binned in order to have at least fifty photons for 
each bin. Obviously, both re-binning are made such that 
minimum chisquare techniques could be reliably used in 
spectral fitting. Spectral analysis was performed using the 
energy range: 1.65— 10.5 keV for the MECS and 15— 60keV 
for the PDS (data in the 60-100 keV PDS range was bad), 
and taking into account the off-axis position of the source. 
We used off-axis response matrices made using the pre- 
cise estimation of EXO 053109-6609.2 off axis. Note that 
by using on-axis spectral response matrices for a source 
that is off-axis, one wrongly estimates the source spec- 
tral parameters. In particular, the flux is underestimated 
while the spectrum appears softer than what it actually 
is (Mcrcghctti ct al. 2000). Correspondingly, any extrap- 
olation towards high-energies may be strongly underesti- 
mated. All energy bins the count rate of which was con- 
sistent with zero were not used in the spectral analysis. 

3. Timing analysis 

We search for coherent pulsations of the source in both 
states. In the low state (second observation), in a period 
interval around the ROSAT period, the power spectrum 
and the epoch folding search show an evident signal at 
v <~ 0.0732 Hz with a high significant peak (~ 10a) corre- 
sponding to a period of ~ 13.67s (see Fig.l). Hereafter, if 
not specified, all errors are la confidence level. To rely 
upon a more accurate period measurement, we used a 
phase fitting technique by fitting the phases of the mod- 
ulation over four intervals of ~ 21000 s each. With this 
method we found a spin period of 13.6705 ± 0.0004 s. 

Timing analysis of the first observation was already 
performed giving a spin period value of 13.67590(8) s 
(Burderi et al. 1998). Applying the same method to 
the March 1997 observation, we found a spin period of 
13.67630 ± 0.00003 s (consistent with the period found 
for the same observation in Burderi et al. 1998) 1 . In the 
latter case, owing to a higher statistic, we looked at the 
phase modulations over six intervals of ~ 14000 s: fitting 
the results with a polynomial function of the second or- 
der we found evidence for a quadratic component pos- 
sibly related to the orbital motion of the X-ray pulsar 
(see Fig. 2). The observation is unfortunately too short 
to look for an orbital solution. From the periods mea- 
sured in these two BeppoSAX observations we inferred a 
P = -(1.14 ± 0.08) x lO^ss" 1 . 

1 Note that while the period in Burderi et al. 1998 is referred 
to the start of the observation, here we report the periods at the 
half time of the observation. Taking in account the local period 
derivative reported before, the two periods are consistent. 
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The pulse profile obtained when the source was in 
high-state (March 1997) is dominated by the first har- 
monic. This predominance disappear at fainter fluxes 
(October 1998) where the fundamental is the strongest 
component (Fig. 3). Furthermore, the pulsed fraction (PF) 
in the MECS range changes from a value of 50 ± 1% to 
30 ± 1% from high to low state, while in the PDS range 
only in the high state we see the pulsation, and the PF 
is 40 ± 4% (see Fig.3). The low state PDS epoch folding 
at the period found in the MECS data did not show any 
modulation due to poor statistics (with an upper limit on 
the PF is <~ 80% or 100% depending if we estimate the 
PDS flux of LMC X-4 be 3.5 counts/s or 4.0 counts/s), 
correspondingly we cannot exclude an high-energy pulsa- 
tion of the source in low state. 

Table 1. Spectral parameters of EXO 0531-66 in high- 
state (absorbed blackbody plus a power-law and a 
Gaussian) and in low-state (absorbed power-law). Fluxes 
are unabsorbed and in unit of xl0~ n ergcm~ 2 s _1 (see 
also Fig. 4); errors are la confidence level. 



Spectral parameters 


High state 


Low state 


^(xlO^cm-^) 


0.92 ± 0.04 


0.92 (frozen) 


Power-law PI 


1.75 ±0.08 


0.4 ±0.1 


flux PL (2-10 keV) 


2.0 




kT bb (keV) 


3.53 ±0.3 




flux bb (2-10 keV) 


1.9 




E Fe (keV) 


6.8 ±0.1 




o Fe (keV) 


0.4 ±0.1 




EQW Fe (eV) 


275 




x, 2 


1.18 (169 dof) 


1.2 (25 dof) 


Total Flux (2-10 keV) 


5.0 


0.5 


Total Flux (2-60 keV) 


14.0 





2 5 10 20 50 

Energy (keV) 




2 5 10 
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Fig. 4. BeppoSAX spectra of EXO 053109-6609.2 in high 
state (1997) and in low state (1998). See Tab.l for spectral 
parameters. 

is improved from 1.28 (172 dof) to 1.18 (169), giving an F 
test probability of 3.5 x 10~ 4 . 

The low-state spectrum is modelled by an absorbed 
power-law (see Tab.l and Fig. 4) and the source luminosity 
in this state is 1.5 x 10 36 crgs _1 (for a distance of 50 kpc; in 
2-60 keV energy range). We also tried to add a Gaussian 
function at ~ 6.8 keV but it was no significative (confi- 
dence level < 1.5cr). 



4. Spectral analysis 

EXO 053109-6609.2 high-state spectrum is well fitted by 
an absorbed blackbody plus power-law component and 
a Gaussian function (see Tab.l and Fig. 4). The inferred 
luminosity between 2-60 keV, assuming a distance of 50 
kpc (this is the distance known to the Large Magellanic 
Cloud), is 4.2 x 10 37 ergs _1 . We tried to fit the spec- 
trum with other models: two blackbody plus a power- 
law, a blackbody plus a cutoff power-law and a black- 
body plus a high energy cutoff model but all gives worst 
results (xt > 1-9)- Actually, fitting the spectrum with a 
blackbody plus a cutoff power-law we initially obtained a 
good chisquare value but the absorption was too small ( 
N H = (0.03 ±0.01) x 10 22 atoms/cm 2 ) for a system in the 
LMC. We then fixed the absorption value at ~ 0.1 x 10 22 
atoms/cm 2 but the chisquare becomes remarkably worse 
(this is the galactic absorption value calculated for the RA 
and DEC of EXO 053109-6609.2 with the nh ftool; Dickey 
& Lockman, 1990). Adding a Gaussian line at ~ 6.8 keV 
(see Tab. 1) to the continuum model the chisquare value 



5. Discussion 

During the outburst phase of Be X-ray binaries, accre- 
tion disks are expected to be present, and indeed, evidence 
for an accretion disk, based on a correlation between the 
observed flux and spin-up or spin-down rate , has been 
found for several sources (Bildstcn ct al. 1997; Okazaki 
& Negueruela 2001). The EXO 053109-6609.2 spin period 
was detected four times: the first time by ROSAT (Dennerl 
et al. 1996), two times by BeppoSAX (Burderi et al. 1998 
and this paper) and the fourth time by XMM-Newton 
satellite (Haberl et al. 2003). In Fig.5 we plot the 0.1- 
2.4 keV band luminosity versus period for all these four ob- 
servations; it is evident from the plot that periods and lu- 
minosities are directly correlated. Comparing the spin pe- 
riods of all the observations of EXO 053109-6609.2 carried 
out from its discovery, we found that the source alternates 
spin-up and spin-down (Fig.5). The secular period deriva- 
tive between the ROSAT observation (Haberl, Dennerl, 
and Pietsch 1995) and the first BeppoSAX observation 
is (2.9 ± 0.1) x 10 _11 ss _1 ; comparing the two BeppoSAX 
observations we obtain — (1.14±0.08) x 10~ 10 ss _1 , and the 
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Time (MJD) 



Fig. 5. Spin period and luminosity changes between ROSAT (Haberl et al. 1996), BeppoSAX (this paper) and XMM 
(Haberl et al. 2003) observations (all reported luminosities are extrapolated in 0. 1-2.4 keV band). 



P between the last BeppoSAX observation and the XMM 
observation (Haberl et al. 2003) is -(3.7±0.1) x lO" 1 ^" 1 
(period derivative errors are at 90% confidence level). The 
secular spin period derivative is changed from a spin-down 
to a spin-up trend, which has a clear correlation with the 
X-ray luminosity. We detected this change around 1997. 

Double peaked pulse profile (Fig. 3, first and second 
panels) was observed when the source was in high-state, 
while a single peaked profile was found in the low-state 
(Fig. 3, third panel). This behaviour has been often seen in 
X-ray binary systems and has been ascribed to the tran- 
sition of the source between a pencil beam and a pencil 
plus fan beam emission geometry, a behaviour which is 
correlated with the X-ray flux (Parmar et al. 1989a and 
1989b). Timing behaviour of this source is similar to what 
seen in other Be X-ray binaries, making the period deriva- 
tive changes a peculiar characteristic of these systems. 

The spectrum in the high-state of emission is well 
fitted by an absorbed blackbody plus a power law and 
a Gaussian line at ~ 6.8 keV (see Tab.l and Fig. 4 first 
panel). The 6.8 keV emission line is probably due to K- 
shcll emission from highly ionized iron (probably in the 
He-like and/or H-like ionization stages). This is the first 
evidence of an iron emission line in this source; iron emis- 
sion lines arc a common feature in HMXBs, although their 
centroid energies are usually detected in the range 6.4 - 
6.7 keV (e.g. Parmar et al. 1989a). 

From measured temperature of kT^ 3.5 keV we in- 
fer an emission radius for the blackbody component of 
~ 5.5 km, about half neutron star radius. The calculated 



blackbody radius is smaller than the neutron star radius, 
probably due to the fact that the thermal emission does 
not come from all neutron star surface but probably just 
from a small region around the polar caps. The energy flux 
of the two spectral components differs of about two orders 
of magnitude, Fbb/F power — 1.728 x 10~ 2 . The blackbody 
component probably disappear while the system is in the 
low-state and the power-law photon index becomes flatter 
(see Tab.l and Fig. 4 second panel). 
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